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1. Introduction 

The c = 1 string tlieoryS lias been originally constructed in the early 90 's as the 
theory of random surfaces embedded into a one- dimensional spacetime 0. Since then it 
became clear that this is only one of the realizations of a more universal structure, which 
reappeared in various mathematical and physical problems. Most recently, it was used 
for the description of the 2D black hole |^,H,|5|J|] in the context of the topological string 
theories on Calabi-Yau manifolds with vanishing cycles [0 and M = 1 SYM theories . 

The c = 1 string theory can be constructed as the collective field theory for a one- 
dimensional N X N hermitian matrix field theory known also as Matrix Quantum Mechan- 
ics (MQM). This construction represents the simplest example of the strings/matrix corre- 
spondence. The collective excitations in the singlet sector of MQM are massless "tachyons" 
with various momenta, while the non-singlet sectors contain also winding modes. 

The singlet sector of MQM can be reduced to a system of nonrelativistic fermions 
in the upside-down gaussian potential. Thus the elementary excitations of the c = 1 
string can be represented as collective excitations of free fermions near the Fermi level. 
The tree-level string-theory 5'-matrix can be extracted by considering the propagation of 
infinitesimal "pulses" along the Fermi sea and their reflection off the "Liouville wall" P,p!0|. 



The (time-dependent) string backgrounds are associated with the possible profiles of the 
Fermi sea . 



An important property of the c = 1 string theory is its integrability. The latter has 
been discovered by Dijkgraaf, Moore and Plesser |T^] in studying the properties of the 
tachyon scattering amplitudes. It was demonstrated in fl^ that the partition function of 



the = 1 string theory in the case when the allowed momenta form a lattice, as in the 
case of the compactified Euclidean theory, is a tau function of Toda hierarchy [|13|. The 
operators associated with the momentum modes in the string theory have been interpreted 
in [jl2| as Toda flows. A special case represents the theory compactified at the self-dual 
radius R = 1. It is equivalent to a topological theory that computes the Euler characteristic 
of the moduli space of Riemann surfaces [|14[. When R = 1 and only in this case, the 



partition function of the string theory has alternative realization as a Kontsevich-type 



model p|,|Tl]. 



In this paper we will show that there is another realization of the c = 1 string theory 
by the so-called Normal Matrix Model (NMM). This is a complex N x N matrix model in 
which the integration measure is restricted to matrices that commute with their hermitian 
conjugates. The normal matrix model has been studied in the recent years in the context of 
the laplacian growth problem, the integrable structure of conformal maps, and the quantum 



Hall droplets |[Tq,pr7|,p!8|jT9|,pO ,21]. It has been noticed that both matrix models, the singlet 



sector of MQM and NMM, have similar properties. Both models can be reduced to systems 



^ As a good review we recommend 
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of nonrelativistic fermions and possess Toda integrable structure. There is however an 
essential difference: the normal matrix model describes a compact droplet of Fermi liquid 
while the Fermi sea of the MQM is non-compact. Furthermore, the perturbations of NMM 
are introduced by a matrix potential while these of MQM are introduced by means of 
time-dependent asymptotic states. In this paper we will show that nevertheless these two 
models are equivalent. 

The exact statement is that the matrix quantum mechanics compactified at radius R 
is equivalent to the normal matrix model defined by the probability distribution function 
exp[^WR{Z, Z% where 

WniZ, Zt) = tv{ZZ^f - h^-j^tv\og{ZZ^) - ^t^trZ" - 5^^_fctrZt^ (u) 

A;>1 A;>1 

More precisely, the grand canonical partition function of MQM is identical to the canonical 
partition function of NMM at the perturbative level, by which we mean that the genus 
expansions of the two free energies coincide. This will be proved in two steps. First, we 
will show that the non-perturbed partition functions are equal. Then we will use the fact 
that both partition functions are r-functions of Toda lattice hierarchy, which implies that 
they also coincide in presence of an arbitrary perturbation. 

We also give a unified geometrical description of the two models in terms of a complex 
curve with the topology of a sphere with two punctures. This complex curve is analog 
of the Riemann surfaces arising in the case of hermitian matrix models. However, for 
generic R it is not an algebraic curve. The curve has two dual non-contractible cycles, 
which determine the boundaries of the supports of the eigenvalue distributions for the 
two models. The normal matrix model is associated with the compact cycle while the 
MQM is associated with the non-compact cycle connecting the two punctures. We will 
construct a globally defined one-form whose integrals along the two cycles give the number 
of eigenvalues and the derivative of the free energy with respect to A^. 

The paper is organized as follows. In the next section we will remind the realization of 
the c = 1 string theory as a fermionic system with chiral perturbations worked out in ||TI 



We will stress on the calculation of the free energy, which will be needed to establish the 
equivalence with the NMM. In Sect. 3 we construct the NMM having the same partition 
function. In Sect. 4 we consider the quasiclassical limit and give a unified geometrical 
description of the two models. 
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2. Matrix quantum mechanics, free fermions and integrability 

2.1. Eigenf unctions and fermionic scattering 

In absence of winding modes, the 2D string theory is described by the singlet sector 
of Matrix Quantum Mechanics in the double scaling limit with Hamiltonian 

Ho = ^tr(-n^^-X^). (2.1) 

The radial part of this Hamiltonian is expressed in terms of the eigenvalues xi, . . . ,xn of the 
matrix X. The wave functions in the S'C/(Ar)-singlet sector are completely antisymmetric 
and thus describe a system of N nonrelativistic fermions in the inverse gaussian potential. 
The dynamics of the fermions is governed by the Hamiltonian 



1 ^ 

Ho-^B^^'-^')' (2-2) 



where pi arc the momenta conjugated to the fermionic coordinates Xi. 

To describe the incoming and outgoing tachyonic states, it is convenient to introduce 
the "light-cone" coordinates in the phase space 



satisfying the canonical commutation relations 

[x^,x_] = -in. (2.4) 

In these variables the one-particle Hamiltonian takes the form 

Ho = ~{x^x_ + x_x^). (2.5) 

We can work either in x_^_ or in x_ representation, where the theory is defined in terms of 
fermionic fields V'i (a^± ) respectively. The solutions with a given energy are {x_^ , t) = 

where the phase factor (j)Q{E) will be determined below. These solutions form two complete 
systems of 5-function normalized orthonormal states under the condition that the domain 



e-i-^*V'+(a^±) with 



il^lix^) = -j= e^^'^°a;^^^-5, (2.6) 
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of the definition of the wave functions is a semi-axisJl The left and right representations 
are related to each other by a unitary operator which is the Fourier transform on the 
half-line 

/"OO 

[SiIj+]{x_)= dx^ S{x^,x^)'iIj_^_(x_^). (2.7) 







The integration kernel S{x_,x^) is either a sine or a cosine depending on the boundary 
conditions at the origin. For definiteness, let us choose the cosine kernel 

1 




S(X_,X:) = \ —COS(—X:X). 

Tvn Ti 



Then the operator S is diagonal on the eigenf unctions ( |2.6| ) with given energy 

[Si)l\{x_) = e-^'^o7^(E)^/>f (a;_), (2.9) 

where 

n{E) = TC^\l^ cosh r (^^E+1/2^ . (2.10) 

Since the operator 5" is unitary, i.e. 



n{E)n{E) = n{-E)n{E) = i, (2.11) 

it can be absorbed into the function (^q{E). This fixes the phase of the eigenfunctions as 

0o(S) = -in\ogn{E). (2.12) 

In fact, the operator S appears in our formalism as the fermionic S'-matrix describing 
the scattering off the inverse oscillator potential and the factor TZ{E) gives the fermionic 
reflection coefficient. 

Let us introduce the scalar product between left and right states as follows 



{4>^\S\ij^) = jj dx_^dx__ ip_{x_)S{x_,x^)ip_^{x^). (2-13) 



^ In order to completely define the theory, we should define the interval where the phase- 
space coordinates are allowed to take their values. One possibility is to allow the eigenvalues 
to take any real value. The whole real axis corresponds to a theory whose Fermi sea consists 
of two disconnected components on both sides of the maximum of the potential. To avoid the 
technicalities related to the tunneling phenomena between the two Fermi seas, we will define the 
theory by restricting the eigenvalues to the positive real axis. The difference between the two 
choices is seen only at the non-perturbative level. 
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Then the in- and out-eigenfunctions are orthonormal with respect to this scalar product 



(:iP''_\S\iPl') =5{E ~ E'). (2.14) 

It is actually this relation that defines the phase factor containing all information about 
the scattering. 

The ground state of the Hamiltonian (|2.5| ) can be constructed from the wave func- 
tions (^]^) and describes the linear dilaton background of string theory [ TH ] . The tachyon 



perturbations can be introduced by changing the asymptotics of the wave functions at 

x± ^ oo to 

*|(xj ~ e^^^^*^''<'\±^^-i (2.15) 

The exact phase contains also a constant mode ^^(f){E) as in ( |2.6|) and negative powers 
of a;^/^ which are determined by the orthonormality condition 

{^^\S\^f) = S{E-E'). (2.16) 

The constant mode (p{E) of the phase of the fermion wave function contains all essential 
information about the perturbed system. 

2.2. Cut- ojf prescription and density of states 

To find the density of states, we introduce a cut-off A by confining the phase space to 
a periodic box 

x^+2VA = x^. (2.17) 



which can be interpreted as putting a refiecting wall at distance vA. This means that at 
the points = and x_ = V7i the refiected wave function coincides with the incoming 
one 

[^^'](VA) = ^'(VA). (2.18) 

Applied to the wave functions ( |2.15| ), this condition gives an equation for the admissible 
energies 

^^HE) ^ g--y(A)^tE^ y^^^ ^ j^J2{tk + t-.)A'=/2« (2.19) 

where we neglected the negative powers of A. It is satisfied by a discrete set of energies 
En defined by the relation 

EnlogA-(p{En) ^2nhn + V{A), n e Z. (2.20) 

Taking the limit A — oo we find the density of states (in units of h) 

2tv 2tv dE 
5 



2. 3. The partition function of the compactified 2D string theory 

Knowing the density of states, we can calculate the grand canonical partition function 
Z of the quantum-mechanical system compactified at Euclidean time /3 = 27rR. If fi is the 
chemical potential, then the free energy = h \ogZ of the ensemble of free fermions is 
given by 

/oo 
dEp{E)\og l + e-*^('^+^) . (2.22) 
- oo 

Integrating by parts and dropping out the A-dependent non-universal piece, we obtain 

We close the contour of integration in the upper half plane and take the integral as a sum 
of residues of the thermal factori 

n>0 ^ ^ 

We will be interested only in the perturbative expansion in l/^u and neglect the non- 
perturbative terms ~ e~^'^, as well as non-universal terms in the free energy (regular in 
/i). Therefore, for the case of zero tachyon couplings we can retain only the F- function in 
the reflection factor ( 2.1U ) and choose the phase (f)Q as 

ME) = -in log r (^-E + 1/2^ . (2.25) 

It is known []T^ that the tachyon scattering data for MQM are generated by a r- 
function of Toda lattice hierarchy where the coupling constants tk play the role of the 
Toda times. The r-function is related to the constant mode (t){E) by |11] 

e^^(-) = !4^i±!Al. (2.26) 

Taking into account ( p.24|) we conclude that the partition function of the perturbed MQM 
is equal to the r-function: 

Z{ii,t) = To{ii,t). (2.27) 

The discrete space parameter s G 2 along the Toda chain is related to the chemical 
potential p. More precisely, s corresponds to an imaginary shift of n P,PT|: 



rs(/u,t) =ro + . (2.28) 

This fact was used to rewrite the Toda equations as difference equations in ii rather than 
in the discrete parameter s. It will be also used below to prove the equivalence of 2D string 
theory to the normal matrix model. 



We neglect the nonperturbative terms associated with the cuts of the function (f>(E). 



6 



3. 2D string theory as normal matrix model 

In this section we will show that the partition function of 2D string theory with 
tachyonic perturbations can be rewritten as a normal matrix model. It is related to the 
original matrix quantum mechanics by a kind of duality transformation. In fact, we will 
give two slightly different normal matrix models fulfilling this goal. 

Consider the following matrix integral 

Zr"{N,t,a)= J dZdZ^ [det(ZZt)]^''-'^/'+*e-^*'^[(^^')"-^+(^)-^-(^')], (3.1) 

[Z,Zi]=0 

where the integral goes over all complex matrices satisfying [Z, Z^ = and the potentials 
are given by 

V±iz) = J2i±kz''. (3.2) 

k>i 

We made the dependence on the Planck constant explicit because later we will need to 
analytically continue in h. The integral (|3.1| ) defines a r-function of Toda hierarchy |l20[] . 



We remind the proof of this statement in Appendix A, where we also derive the string 
equation specifying the unique solution of Toda equations. This string equation coincides 
(up to change h —>■ i%) with that of the hierarchy describing the 2D string theory |]6|,|TT 



Therefore, if we identify correctly the parameters of the two models, the r-functions should 
also be the same. Namely, we should relate the chemical potential /i of 2D string theory 
to the size of matrices and the parameter a of the normal matrix model. There are two 
possibilities to make such identification. 

3.1. Model I 

The first possibility is realized taking a large N limit of the matrix integral ( |3.1|) . 
Namely, we will prove that the full perturbed partition function of MQM is given by the 
large N limit of the partition function ( |3.1| ) with a = R^i — %N and a subsequent analytical 
continuation % ^ i% 

Zn{ix,t)= lim zX^{N,t,Rfi-inN). (3.3) 

A''— >oo 

The necessity to change the Planck constant by the imaginary one follows from the com- 
parison of the string equations of two models as was discussed above. More generally, 
the T-function (|2.28 ) for arbitrary s is obtained from the partition function ( |3.1| ) in the 
following way 



r,,n(^,i)= lira ZZ'^iN + s^t^R^-ihN). (3.4) 

Let us stress that in spite of the large N limit taken here, we obtain as a result the full 
(and not only dispersionless) partition function of the c = 1 string theory. 
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To prove this statement we will use the fact that both partition functions are r- 
functions of the Toda lattice hierarchy with times t±k, /c = 1, 2, . . .. Since the unperturbed 
r-function provides the necessary boundary conditions for the unique solution of Toda 
equations, it is sufficient to show that the unperturbed partition functions of the two 
models coincide and the s-parameters of two r-functions are identical. 

The integral reduces to the product of the normalization coefficients of the or- 
thogonal polynomials, eq. (A. 5). From the definition of the scalar product (A. 4) it follows 
that when all t±k = the orthogonal polynomials are simple monomials and the normal- 
ization factors are given by 

h^(^O^a) = — [ d^ze-^^''^''{zzY''-^^M+^. (3.5) 
2m Jc 

Up to constant factors and powers of we find 

We see that the result coincides with the reflection coefficient given in (|2.25| ). Therefore, 
taking the analytical continuation of parameters as in (|3.3| ), we have 

^-1 / AT ^ 1 

loe ZT^ {N, 0, Ru - ihN) = lim ih 

AT— >oo 



lim logZT'' (iV, 0, Ru - ihN) = hm V (/)o ih ^ - n 

n=0 ^ ^ 

°° f n + - \ 



This coincides with the unperturbed free energy J-'{n, 0) from ( p.24|) , what proves (|3.3|) for 
all t±k = 0. 

Thus, it remains to show that the s-parameter of the r-function describing the normal 
matrix model is associated with fj,. From (A. 17) we see that it coincides with A^. Hence it 
trivially follows from (|3.4D that 



Ts,n{f^, t) = lirn^ ZT"^ (a, t, i? + - z^a) = To,n [fi + ih^, t) , (3.8) 

what means that the r-function defined in (^]J) possesses the characteristic property ( p.28|) . 



3.2. Model II 

In fact, one can even simplify the representation ( |3.3| ) what will provide us with the 
second realization of 2D string theory in terms of NMM. Let us consider the matrix model 
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for a = 0. We claim that the partition function Zifi{N, t, 0) in the canonical ensemble, 
analytically continued to imaginary Planck constant and 

N = -^Rfi, (3.9) 
n 

coincides with the partition function of 2D string theory in the grand canonical ensemblefl 

ZnM=z:^i~R^^.t). (3.10) 

Indeed, as for the first model, they are both given by r-functions of Toda hierarchy. 
After the identification ( p.9| ), the s-parameters of these r-functions are identical due to 
(A. 17) and ( p.28|) . Therefore, it only remains to show that the integral Zifi{N, 0) without 
potential is equal to the unperturbed partition function ( p.24| ). In this case the method of 
orthogonal polynomials gives 

N-l y , 1 



zT^{N,o)=llnUn^^), (3.11) 

n=0 ^ ^ 

where the TZ factors are the same as in MQM, as was shown in (^). Then we write 

^NMM ^ S(0)/S(Ar), where E{N) = f] 7^ f-zHn + ^)/i?) . (3.12) 

n=N ^ ^ 

S(0) is a constant and can be neglected, whereas S(A^) can be rewritten as 

E{N) = n 7^ (-ihN/R - ih{n + -) /r) . (3.13) 

n=o ^ 2 y 

Taking into account the unitarity relation (|2.11| ) and substituting from ( p.9|) , we obtain 

^T^ (-^^^, 0) ~ s-H-^i?^) = f[nL + in{n + \)/r) . (3.i4) 

n=0 ^ ^ 



Thus fi^ log Z^j^'^ { — ^R^,0) coincides with the free energy ( p.24|) with all tk = 0. Note 
that the difference in the sign of ^ does not matter since the free energy is an even function 
of n (up non- universal terms) .H Since the two partition functions are both solutions of the 



Prom now on we will omit the last argument of Z^^^ corresponding to the vanishing pa- 
rameter Q. 

^ Moreover, this sign can be correctly reproduced from ( p.23| ) if closing the integration contour 
in the lower half plane. 
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To da hierarchy, the fact that they coincide attk = imphes that they coincide for arbitrary 
perturbation. 

We conclude that the grand canonical partition function of the c = 1 string theory 
equals the canonical partition function of the normal matrix model ( |3.1| ) for a = and 
= —^Rfi. It is also clear that we can identify the operators of tachyons in two models 

TrXl^^ ^ TrZ"", Tr ^ ^ Tr {Z^ . (3.15) 



4. Geometrical meaning of the duality between the two models 

The quasiclassical limit of most of the solvable matrix models has a nice geometrical 
interpretation. Namely, the free energy in this limit can be parameterized in terms of the 
periods of a holomorphic 1-form around the cycles of an analytical curve [ p2| , p3| , |2^ , p5 1 . 



Recently this geometrical picture appeared also in the context of the topological strings 
on singular Calabi-Yau manifolds and supersymmetric gauge theories [I26|j2^ , p8| . 

In this section we will show that both MQM and NMM have in the quasiclassical limit 
a similar geometrical interpretation in terms of a one-dimensional complex curve, which 
is topologically a sphere with two punctures. Here by complex curve we understand a 
complex manifold with punctures and given behavior of the functions on this manifold at 
the punctures. Each of the two matrix models corresponds to a particular real section of 
this curve coinciding with one of its two non-contractible cycles, and the duality between 
them is realized by the exchange of the cycles of the curve. 

4.1. The dispersionless limit 

The quasiclassical limit h corresponds to the dispersionless limit of the To da 
hierarchy where it has a description in terms of a classical dynamical system. The Lax 
operators are considered as c-functions of the two canonically conjugated coordinates fi 
and where u is the quasiclassical limit of the shift operator uj = e^l^^ . Moreover, the 
two Lax operators can be expressed as series in u). The solutions of the dispersionless 
hierarchy correspond to canonical transformations in the phase space of the dynamical 
system. In the case of MQM this is the transformation from the coordinates \i and log uj 
to and x„ . 

The particular solution of the Toda hierarchy that appears in our problem also satisfies 
the dispersionless string equation. In the case of MQM this is nothing but the equation of 
the profile of the Fermi sea in the phase space [0 

x^x_ = ^ kt±kX^J^ + + 5Z v±kxl^'^- (4-1) 
fc>i fc>i 

To get the string equation for NMM it is enough to make the substitution following ( |3.15| ) 

x^ ^ z^, x^ ^ (4.2) 

and n = hN/R as explained in the previous section. (We took into account that h from 
( P^ ) should be replaced here by —ih so that the factor i is canceled.) Then the string 
equation describes the contour 7 bounding the region T> in the complex ^-plane filled by 
the eigenvalues of the normal matrix.! 



Here we consider the case when the eigenvalues are distributed in a simply connected domain. 
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4-. 2. NMM in terms of electrostatic potential 

In this section we will introduce a function of the spectral variables z and z, which plays 
a central role in the NMM integrable structure. This function has several interpretations in 
the quasiclassical limit. First, it can be viewed as the generating function for the canonical 
transformation mentioned above, which maps variables /i and log uj to the variables z and 
z. Second, it gives the phase of the fermion wave function at the Fermi level after the 
identification (^7^). 

There is also a third, electrostatic interpretation, which is geometrically the most 
explicit and which we will follow in this section.0 According to this interpretation, the 
eigenvalues distributed in the domain V can be considered to form a charged liquid with 
the density 



1 „ 



p{z, z) = -d.d.Wniz, z) = — (4.3) 

TT TT 

Let ip{z, z) be the potential of the charged eigenvalue liquid, which is a harmonic function 
outside the domain V and it is a solution of the Laplace equation inside the domain 

To fix completely the potential, we should also impose the asymptotics of the potential 
at infinity. This asymptotics is determined by the coupling constants n = 1, 2, ... and 
can be considered as the result of placing a dipole, quadrupole etc. charges at infinity. 

The solution of this electrostatic problem is obtained as follows. The continuity of the 
potential (fi{z, z) and its first derivatives leads to the following conditions to be satisfied 
on the boundary 7 = dV 

ip{z)+^{z) = {zz)'', (4.5) 
zd,^{z) = zd^^iz) = Rz^z^. (4.6) 
Each of two equations p. 6]) can be interpreted as an equation for the contour 7. Since we 



obtain two equations for one curve, they should be compatible. This imposes a restriction 
on the holomorphic functions fiz) and <f>{z). The solutions for these chiral parts of the 
potential can be found comparing (|4.6| ) with the string equation (|4.1|) . In this way we have 



^{z) = hN log z + ^(j) + Rj^ikz'' -Rj^rVkZ 



k 

Z 

(4.7) 



^{z) = hN log z + ^(I) + rJ2 t-kz'" -RYI \^-^^ 



k 

k>i k>i 



^ See also [29| for another interesting interpretation which appears useful in string field theory. 



11 



The zero mode (j) is fixed by the condition (O). However, it is easier to use another 
interpretation. As we mentioned above, the holomorphic functions f{z) and v^(^) coincide 
with the phases of the fermionic wave functions and in particular their constant modes are 



the same. Then, taking the hmit ^ ^ of ( 2.26 ), in the variables of NMM we find 



-nAlogZ— . (4.^ 



4-3. The complex curve 



Instead of considering the function (fi{z, z) harmonic in the exterior domain f> = C/P, 
we can introduce one holomorphic function $ defined in the double cover of T>. The doubly 
covered domain V is topologically equivalent to a two-sphere with two punctures (at 2 = oo 
and z = 00), which we identify with the north and the south poles. It can be covered by two 
coordinate patches associated with the north and south hemispheres and parameterized 
by z and z. They induce a natural complex structure so that the sphere can be viewed 
as a complex manifold. The patches overlap in a ring containing the contour 7 and the 
transition function between them z{z) (or z{z)) is defined through the string equation 

(E§)- 

The string equation can also be considered as a defining equation for the manifold. 
We start with the complex plane with fiat coordinates z and z. Then the manifold is 
defined as solution of equation (|4.6|) , thus providing its natural embedding in C^. 

To completely define the complex curve, we should also specify the singular behavior 
at the punctures of the analytic functions on this curve. We fix it to be the same as the 
one of the holomorphic parts of the potential ( [4.7| ). 

We define on the complex curve the following holomorphic field 

^def f ^+{z) = Lp{z) — ^{zz{z))^ in the north hemisphere, ^ 
1 ^-{z) = —(fiiz) + ^{z{z)z)^ in the south hemisphere. 



Its analyticity follows from equation ( [4.5| ), which now holds on the entire curve, since z 
and z are no more considered as conjugated to each other. 

The field $ gives rise to a closed (but not exact) holomorphic 1-form globally 
defined on the complex curve. Actually, this is the unique globally defined 1-form with 
given singular behavior at the two punctures. Therefore, we can characterize the curve 
by periods of this 1-form around two conjugated cycles A and B defined as follows. The 
cycle A goes along the ring where both parameterizations overlap and is homotopic to the 
contour 7. The cycle S is a path going from the puncture at 2; = 00 to the puncture at 
z = 00. 



The integral around the cycle A is easy to calculate using (|4.7|) and ( ^1.6|) by picking 
up the pole 

^ ^ d<^ = — (f d<^+{z) = — / d<^-{z) = hN, (4.10) 
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where we indicated that in ^-coordinates the contour should be reversed. 

To find the integral along the non-compact cycle B, one should introduce a regular- 
ization by cutting it at z = y/X and z = Then we obtain 

d^= d^+{z)+ / = $_(yA) -$+(v^), (4.11) 

Jb Jy/X Jz(zo) 

where zq is any point on the cycle and we used ( [4.5| ). Taking into account the definition 
( ^^ ) and the explicit solution ( [4.7| ), we see that the result contains the part vanishing 
at A ^ oo, the diverging part which does not depend on and, as non-universal, can 
be neglected, and the contribution of the zero mode (p. Since (p enters $_|_ and $_ with 
different signs, strictly speaking, it is not a zero mode for $. As a result, it does not 
disappear from the integral, but is doubled. Finally, using ( [1.8|) we get 



^ci$ = nAiogZ^--. (4.12) 
In Appendix B we derive this formula using the eigenvalue transfer procedure similar to 



one used in [^, or in in a similar case of the two matrix model equivalent to our 
R = 1 NMM. Thus, the free energy of NMM can be reproduced from the monodromy of 
the holomorphic 1-form around the non-compact cycle on the punctured sphere. 

Since the string equations coincide, it is clear that the solution of MQM is described 
by the same analytical curve and holomorphic 1-form. To write the period integrals in 
terms of MQM variables, it is enough to change the coordinates according to (|4.2| ) and hN 
= Rfx in all formulas. In this way we findi 

^(f d^ = R^, [ d^ = -i^. (4.13) 
2m J A Jb R oil 

The last integral can be also obtained from the integral over the Fermi sea. Indeed, we 
can choose a point on the contour of the Fermi sea and split the integral, which actually 
calculates the area of the sea, into three parts 

1 dJ' f /"V^A /--v/A 

— -— = — / dx^dx_= I x_{x^)dx^-[- I x^{x_)dx_ -\- xqx_{xq). (4.14) 

R oil Jf.s. Jxo Jx_{xo) 

Since d+^p = x_{x_^_) and d-(fi = x_^_{x_) (see ([4.6|)), the last expression is equal to the 
integral ( |4.11j ) of (i$ around the (reversed) cycle B. This derivation gives an independent 
check that the free energies of NMM and MQM do coincide. 



To get the second integral from ( [4.12 ), one should take into account that also h ih. 



13 



Fig. 1: Symbolic representation of the complex curve and the two real 
sections along the cycles A and B. The filled regions symbolize the Fermi 
seas of the two matrix models. 

4 ■4- Duality 

We found that the solutions of both models are described by the same complex curve. 
The curve is characterized by a pair of conjugated cycles: the compact cycle A encircling 
one of the punctures and the non-compact cycle B connecting the two punctures. The 
parameter of the free energy fj, or N and the derivative of the free energy itself are given 
by the integrals of the unique globally defined holomorphic 1-form along the A and B 
cycles, correspondingly. 

Now one can ask: what is the difference between the two models? Is it seen at the level 
of the curve? In fact, both models, NMM and MQM can be associated with two different 
real sections of the complex curve. Indeed, in NMM the variables are conjugated to each 
other, whereas in MQM they are real. Therefore, let us take the interpretation where the 
curve is embedded into and consider its intersection with two planes. The first plane 
is defined by the condition z* = z and the second one is given by z* = z, z* = z. In the 
former case we get the cycle A, whereas in the latter case the intersection coincides with 
the cycle B (see fig. 1). One can think about the planes as the place where the eigenvalues 
of NMM and MQM, correspondingly, live (with the density given by ( [4.3| )). Then the 
intersections describe the contours of the regions filled by the eigenvalues. We see that for 
NMM it is given by the compact cycle A, and for MQM the Fermi sea is bounded by the 
cycle B of the curve. 

Therefore, the duality between the two models can be interpreted as the duality 
exchanging the cycles of the complex curve describing the solution. Under this exchange 
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the real variables of one model go to the complex conjugated variables of another model, 
and the grand canonical free energy is replaced by the canonical one. 



This duality can be seen even more explicitly if one rewrites the relations (|4.13[ ) in 
terms of the canonical free energy of MQM defined as F = J-'+hRfiM, where M = 
is the number of eigenvalues. Then they take the following form 

' d^ = -T^, d^ = nM. (4.15) 

Jb 



2ni Ta ndW' 



They have exactly the same form as (|4.10D and ( [1.1 2|) provided we change the cycles. In 



particular, the numbers of eigenvalues are given in the two models as integrals around the 
dual cycles. 



5. Conclusions and problems 

We have shown that the the tachyonic sector of the compactified c = 1 string theory 
is described in terms of a normal matrix model. This matrix model is related by duality 
to the traditional representation in terms of Matrix Quantum Mechanics. More precisely, 
the canonical partition function of NMM is equal to the grand canonical partition function 
of MQM, provided we identify the number of eigenvalues of NMM with the chemical 
potential /U in MQM. This holds for any compactification radius. Moreover, the duality 
between the two matrix models has a nice geometrical interpretation. In the quasiclassical 
limit the solutions of the two models can be described in terms of a complex curve with 
the topology of a sphere with two punctures. The Matrix Quantum Mechanics and the 
Normal Matrix Model are associated with two real sections of this curve which are given by 
the two non-contractible cycles on the punctured sphere. The duality acts by exchanging 
the two cycles. 

The duality we have observed relates a fermionic problem with non-compact Fermi sea 
and continuous spectrum to a problem with compact Fermi sea and discrete spectrum. The 
solutions of the two problems are obtained from each other by analytic continuation. Thus, 
the reason for the duality can be seen in the common analytic structure of both problems. 
We also emphasize that the duality relates the partition functions in the canonical and 
the grand canonical ensembles, which are related by Legendre transform. This is natural 
regarding that it exchanges the cycles of the complex curve. 

Let us mention two among the many unsolved problems related to the duality between 
these two models. First, it would be interesting to generalize the above analysis to the case 
when the eigenvalues of NMM form two or more disconnected droplets. For R = 1 this is 



the analog of the multicut solutions of the two- matrix model, considered recently in |25 



In terms of MQM this situation would mean the appearance of new, compact components 
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of the Fermi sea. The corresponding complex curve has the topology of a sphere with two 
punctures and a number of handles. 

Second, we would like to generalize the correspondence between the two matrix models 
in such a way that it incorporates also the winding modes. For this we should understand 
how to introduce the winding modes in the Normal Matrix Model. Whereas in MQM they 
appear when we relax the projection to the singlet sector, in NMM this could happen when 
we relax the normality condition [Z, Z] = 0. 

Up to now there were two suggestions how to include both the tachyon and wind- 



ing modes within a single matrix model. In [^, a 3-matrix model was proposed with 
interacting two hermitian matrices and one unitary matrix. This model can be seen as a 
particular reduction of Euclidean compactified MQM. The model correctly describes the 
cases of only tachyon or only winding perturbations, though its validity in the general case 
is still to be proven. For the particular case of the self-dual radius and the multiples of it. 



a 4-matrix model was recently proposed in 0. It is based on the old observation of [^| 
about the geometry of the ground ring of c = 1 string theory. It would be interesting to 
find the relation of this model to our approach, at least in the above-mentioned particular 
cases. The general understanding of this important problem is still missing. 
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the RTN contracts HPRN-CT-2000-00122 and -00131. The work of S.A. and I.K was sup- 
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Appendix A. Toda description of the perturbed system 



Here we show that the partition function (|3.1|) is a r-function of the Toda lattice 
hierarchy (for details see We will do it following the standard method of orthogonal 

polynomials. First, we write the integral ( pT|) as integral over the eigenvalues zi, zn 



1 

m 



N 



Y\ WaiZk, Zk)A{z)A{z) 



(A.l) 



k=l 



with measure 



w„iz, z) 



fz_ 
2m 



-i[(zz)«-V+(z)-\/„(5)]^^-^(i?-l)/2+f _ 



(A.2) 
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Then we introduce a set of bi-orthogonal polynomials 

^+(-)=^ / n ?^"^';;';| A(z)A(.-)n(z-z.), (a.3) 



k=l 



and similarly for P^{z). Here A{z) = Yl i^k ~ zi) is the Vandermond determinant and 

k<i 

normalization factors hk are determined by the orthonormality condition 

{P!^\P^).^ = / wUz, z) P!^{z)Pi;j{z) = 5n-m. (A.4) 
Jc 

Then the partition function ( |A.1| ) reduces to the product of the normalization factors 

N-l 

Zr"{NXa)= n /in(i,«). (A.5) 

n=0 

The operators of multiplication by z and z are represented in the basis of orthogonal 
polynomials by the infinite matrices 

ZPI{Z) =J2LnmP^{z), zPI^iz) = Y,P-iz)Lmn. (A.6) 

m m 

with 

Ln,n+1 = -^n+l,n = a/ hn/hn+1, 

(A.7) 

Lnm = Lmn = 0, 771 > U + 1. 

Differentiation of the orthogonality relation ( |A.4| ) with respect to coupling constants tk 
gives 

dP'^(z) """^ 1 

= - J2(L%mP^{z) - -{L^)n.Pl{z), 

dP'^(z) - 1 

m=0 

and similarly for P'^{z). Let us define now a wave function which is a vector \1/ = {\E'n} 
with elements 

^r.{t-z)=Pl{z)e'n^+^'\ (A.9) 
Then the equations ( [A .81 ) lead to the following eigenvalue problem 

z^ = L^, n— = Hk^, n-— = H_k'^, (A.IO) 

Otk Ot_k 

where we introduced the Hamiltonians 

Hk = {L')+ + 1{L%, H_, = -{L^)_ - \{L%. (A.ll) 
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Here the subscripts 0/ — /+ denote diagonal/lower /upper triangular parts of the matrix. 
Prom the commutativity of the second derivatives, it is easy to find the Zakharov-Shabat 
zero-curvature condition 

n^-n^-^ + [H,,H,] = Q. (A.12) 
oil dtk 

It means that the Hamiltonians generate commuting flows and the perturbed system is 
described by the Toda Lattice hierarchy. 

In addition, one can obtain the string equation for the hierarchy. It follows from the 
Ward identity 



z— = -L^L^^f 

oz % 



{R+l)/2+^' ^ (A.13) 
n 



and can be written as 

[L^,L^]=n. (A.14) 

Here we should understand the operators as analytical continuation of the operators 
in integer powers. 

The Toda structure leads to an infinite set of PDE's for the coefficients of the operators 
L and L. The first of these equations can be written for the normalization factors and is 
known as Toda equation 

^ ^^^log^n = 7 7—- (A.15) 

The quantity is known to be related to the r-function of the Toda hierarchy by 

NMM 

hn = (A.16) 
Taking into account ( [A.5|) , this gives for the partition function 

^NMM / AT , \ NMM / , \ / A 1 '7^ 

Zf^ {N,t,a) =r^. {t,a). (A.17) 



Appendix B. Derivation of the formula for the cycle B 

The formula ( ^4.12| ) can be obtained by means of the procedure of transfer of an 
eigenvalue from the point (zi, zi = ^(^i)) belonging to the boundary 7 of the spot to 00 
(or rather to the cut-off point "\/A), worked out for the case R = 1 m p5|]. On the double 



cover described in subsection 4.3, zi and zi are considered as complex coordinates of the 
point in two patches. We find from ( |A.1| ) that the free energy in the large N limit changes 
during this transfer as follows: 



d 



n— logz""" = -(^i^i)^ + + y-(^i) + hY. i°s[(^i - - ^-)]' (B-i^ 
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where all taken at their saddle point values. Here we took into account that the 

determinant in the matrix integral ( |3.1| ) does not contribute in the quasi classical limit. In 
this limit the integral ( |A.1|) leads to the following saddle point equations 



VL{z) = G{z), 



V:,{z) = G{z), 



(B.2) 



where G{z) and G{z) are the resolvents of the one dimensional distributions of z^'s and 
^fc's, respectively .i They have the following asymptotics at large z and z: 



G{z) 



UN 



G{z) 



UN 



(B.3) 



Rewriting the sum in ([B.l|) as an integral with the measure given by the eigenvalue density 
and then expressing it through the resolvents, we find 



n 



d 

dN 



logZ™" = -(zi^i)''+^+(^i)+^-(^i)+/7^G(^)log(zi-^)+/;^G(^)log(zi-z) 

J ZTTl J ZTXI 

(B.4) 

where the contours of integration encircle the whole support of the distribution of eigen- 
values on the physical sheets of the functions z{z) and z{z), respectively.0 

Due to (p.2| ) we do not have any singularities at z oo (z ^ oo) except poles. 
Blowing up the contour we pick up the logarithmic cut and obtain 



^ 1 -7^ 



,ZlZi 



R 



dz 



{zz{z)) 



R 



R 



Zl 



Va 



dz 



f{z{-z)z)^. 



:b.5) 



The first term in the r.h.s. can be regrouped with the other two terms, giving (up to a 
nonuniversal term ~ A-^) 



n 



d 



loffZ 



Zl 



R{zz{z)) 



R 



dz 1 



d{zz 



Zl 



R{z{z)z 



,R 



dz 1 



d{zz 



7^\R 



(B.6) 

which immediately yields if we take into account (j^^ and the definition ( [4.9|) . Note 

that ( [B.2| ) is essentially the same as ( [4.6| ) if one uses the explicit expansion ( |4.7| ) for the 
holomorphic parts of the potential. Therfore, they define the same functions z{z) and z{z). 
On the complex curve described in the subsection 4.3, the formula (|B.6|) reduces to the 
period integral (|4.12| ) of the holomorphic differential d^. 



We view z{z) and z{z) as analytical functions similarly to the interpretation of subsection 4.3. 
Note that the point of infinite branching at the origin related to the singularities of z^ and 
z^ should not appear on the physical sheets. 
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